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STRUCTURAL STRESS ANALYSIS 



CROSS-REFERENCE TO RELATED APPLICATIONS 
The present application claims the benefit of U.S. Provisional Application Serial 
No. 60/249,800 (Attorney Docket No. BAT 0031 MA) for STRUCTURAL STRESS 
CALCULATION SCHEME, filed November 17, 2000. 

BACKGROUND OF THE INVENTION 
The present invention relates to structural stress analysis and, more particularly 
to stress calculations and a scheme for calculating structural stress where the structure 
contains geometric discontinuities, e.g., welded joints, notches, ridges, bends, sharp 
corners, etc. 

Stress analysis is used primarily in the design of solid structures such as ships, 
automobiles, aircraft, buildings, bridges, and dams for determining structural behavior 
and for ascertaining structural viability and integrity under anticipated or foreseeable 
loads. The analysis may involve the use of an abstract or mathematical model for the 
representation of the structural system and loads. According to classical analytical 
idealization, partial differential equations are used. For example, stress in a dam under 
a hydrostatic load can be described by an elliptic partial differential equation in two 
spatial dimensions. 

As boundary geometry of structural systems is usually complicated, the partial 
differential equations of structural mechanics typically cannot be solved in the closed 
analytical exact form. Numerical approximations are sought instead. In one approach, 
derivatives are replaced with finite differences. Other methods are based on finding an 
approximation as a linear combination of preassigned functions such as polynomials or 
trigonometric functions. Also, after a domain or a boundary of interest has been 
discretized in the form of a large number of small elements, a piece-wise approximation 
can be sought according to the finite element method. 

Current methods of stress analysis based upon numeric approximations and 
extrapolation are often subject to substantial uncertainties in regions of close proximity 
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to welds, joints, sudden changes of geometry, or other structural or geometric 
discontinuities and are highly dependent on element size and are typically mesh 
dependent, particularly if drastically different loading modes are considered. 
Accordingly, there is a need for an improved structural stress analysis scheme for 
fatigue prediction that effectively eliminates or minimizes mesh dependency. 

BRIEF SUMMARY OF THE INVENTION 
This need is met by the present invention wherein structural stress in a fatigue- 
prone region of a structure is determined and analyzed by using: 

i) the nodal forces and displacement values in the fatigue- 
prone region, or 

ii) equilibrium equivalent simple stress states consistent with 
elementary structural mechanics in the fatigue-prone region. 

Of course, it is contemplated that combinations, equivalents, or variations of the recited 
bases may alternatively be employed. The determination is substantially independent 
of mesh size and is particularly well-suited for applications where S-N curves are used 
in weld fatigue design and evaluation, where S represents nominal stress or stress 
range and W represents the number of cycles to failure. 

Throughout the present specification and claims certain quantities, functions, 
parameters, and values are described or identified as being determined, calculated, or 
measured. For the purposes of defining and describing the present invention, it is 
noted that the use of one of these terms herein is not intended to exclude steps that 
incorporate or may be defined with reference to another of the terms. For example, a 
"determination" may incorporate aspects that are measured or calculated, a 
"calculation" may depend upon a specific measurement or determination, and a 
measurement may be based upon a calculation or determination. The term "analysis" 
is utilized herein to refer to an identification, review, or examination of the results of a 
determination, calculation, measurment, or combination thereof. 
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For the purposes of defining and describing the present invention, it is noted that 
reference to a "mesh" or "grid" herein relates to the use of any defined reference 
framework for identifying, quantifying, or otherwise defining individual portions of an 
object or an actual or virtual model of an object. Typically, a single quantitative value is 
5 associated with individual points along the reference framework or with individual 

portions defined by the reference framework. A calculation that is mesh insensitive is 
substantially unaffected by the inherent precision, frequency, or definition of the 
reference framework. In stress/strain evaluation, it is common to utilize a finite element 
mesh derived by dividing a physical domain into smaller sub-domains or elements in 
10 J| order to facilitate a numerical solution of a partial differential equation. Surface 

^ domains may be subdivided into triangle or quadrilateral shapes, while volumes may be 

fy subdivided primarily into tetrahedral or hexahedral shapes. The size, shape and 

ill 

£rj distribution of the elements is ideally defined by automatic meshing algorithms. 
^ The present invention is directed to a structural stress calculation scheme that 

1 5 h generates stress determinations that are substantially independent of the size, shape 
and distribution of the mesh elements. Validations of the stress analysis scheme of the 
present invention reveal that structural stress determined according to the present 
invention proved to be substantially independent of mesh size over at least the following 
mesh size ranges: (i) from about 0.1 6* and 0.1/, along respective x and j> axes, to about 
20 2/ and / along respective x and y axes; (ii) from about 0.5/ and /, along respective x and 
y axes, to about 2/ and / along respective x and y axes; and (iii) from about 0.008/ and 
0.02/, along respective x and y axes, to about 0.4/ and 0.5/ along respective x and y 
axes, where / represents the thickness of the structure. For the purposes of defining the 
present invention, it is noted that claims reciting a mesh are intended to cover any 
25 reference framework, like a mesh, grid, or other framework, where a domain is divided 
into sub-domains or elements. 

The present invention has application to stress evaluation of a variety of welded 
and non-welded joints, notches, ridges, bends, sharp corners, and other discontinuities 
or sudden changes in geometry in metallic, plastic and ceramic structures. Relevant 
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structures include aircraft and aerospace equipment (aircraft, helicopters, rockets, 
gantries); agricultural equipment (aerators, balers, baggers, choppers, combines, 
cultivators, elevators, feeders, grain hoppers, bins and sorters, harrows, harvesters, 
mowers, mulchers, planters, plows, scrapers, seeders, shredders, sprayers, spreaders, 
5 tillers, threshers, tractors); agricultural structures (silos, barns, brooders, incubators); 
automobiles and trucks (automobiles, trucks, trailers, wagons, axles, hitches); 
construction and lifting equipment (bulldozers, cranes, hoists, winches, jacks, chains, 
spreaders, hi-los, backhoes, forklifts, loaders, haulers, scrapers, excavators, graders, 
trenchers, borers, directional drillers, pulverizers, saws); forestry equipment (skidders, 
10 J feller bunchers, log loaders, log splitters); mining equipment; rail car frames; ships; 
^ submarines and submersibles; ports and port equipment (docks, piers, cranes, 
PJ lighthouses); bridges (bridges, bridge expansion joints); channels or canals; tunnels; 

in 

^ building components (doors, door frames, windows, sashes, shutters, soffits, fascia, 
^ skylights, scaffolding, cabinetry, chimneys and flues, elevators); building materials 
15 m, (framing channels and struts, joists, trusses and truss plates, wire reinforced concrete); 
appliances, home and industrial (sinks, stoves, ovens, refrigerators, baths); 
5 buildings/skyscrapers (steel decks, mezzanines); housing, particularly manufactured or 
iu modular; heating and cooling systems (especially ducts); home improvement equipment 
(ladders, hand tools); fencing and gates; plumbing (pipes, fittings, pumps, sewage 
20 lines); irrigation and drainage equipment (pipes, irrigation sprinkler systems, drains and 
drainpipes); manufacturing equipment and machinery (conveyors, fasteners, riveters); 
diving equipment; nuclear containers and facilities; offshore oil rigs; diesel and gas 
turbines; pipelines; derricks and digger derricks; cooling towers; radio 
towers/transmitters; welded structures (welded wire reinforcements); tanks and cisterns 
25 (esp. for water storage); aircraft components; automotive parts; footwear; household 
components (sinks, showers, plumbing pipe, swimming pools); sporting goods; 
ceramics; concrete; porcelain enamel; sealants and sealed structures, adhesively 
bonded structures; etc. 

The stress analysis scheme of the present invention may also be utilized to 
30 monitor stress, predict failure, schedule maintenance, determine fatigue, and analyze 
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structures. The structural stress analysis scheme of the present invention is particularly 
significant in these contexts because a 10% variation in a stress value could translate 
into a life cycle that varies by as much as 100-200%. Further, the stress analysis 
scheme of the present invention may be utilized as a structural design tool directed at 
optimizing manufacturing costs or as a tool for monitoring, locally or remotely, structural 
stress in situ or during manufacture or assembly to guide manufacturing or assembly 
steps. Structures may be analyzed by imaging the structure, creating a model of the 
structure and applying the principles of the present invention to assess the stress in the 
structure. It is contemplated that the scheme of the present invention may be 
incorporated into an otherwise conventional, field-ready hand-held or desk top 
computer or programmable device. 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

The following detailed description of the preferred embodiments of the present 
invention can be best understood when read in conjunction with the following drawings, 
where like structure is indicated with like reference numerals and in which: 

Figs. 1A, 1B, 1C, and 1D illustrate typical welded joints and corresponding 
structural stress definitions and local stress distributions associated therewith; 

Figs. 2A and 2B illustrate structural stress definitions and local stress 
distributions associated with spot and laser lap joints; 

Fig. 3 illustrates normal and shear stress distributions from a typical finite 
element model at a reference cross section of a structure for structural stress 
calculations; 

Fig. 4 represents a structural stress definition, taken along section line A-A, 
corresponding to the normal and shear stress distributions of Fig. 3; 

Fig. 5A illustrates local normal stress and shear stress distributions at a 
reference cross-section for a fatigue failure criterion corresponding to a structure 
including partial thickness fatigue crack; 

Fig. 5B represents the structural stress definition relative to the overall thickness 
of the structure of Fig. 5A; 

Fig. 5C represents the structural stress definition for the structure of Fig. 5A with 
a partial thickness fatigue crack failure; 

Fig. 6A illustrates local normal and shear stress distributions of a structure 
defining a non-monotonic through thickness stress distribution at a weld toe of the 
structure; 

Fig. 6B illustrates respective local normal and shear stress distributions at a 
reference section of the structure illustrated in Fig. 6A; 

Fig. 6C illustrates structural stress components relative to both upper and lower 
portions of the structure illustrated in Fig. 6A; 

Fig. 7 is a flow chart illustrating a method of calculating structural stress for 
general three dimensional structures according to the present invention; 
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Fig. 8 illustrates a shell model and corresponding elements for structural stress 
extraction according to one embodiment of the present invention; 

Fig. 9 illustrates a solid model and corresponding elements for structural stress 
extraction according to one embodiment of the present invention; 

Fig. 10 illustrates a scheme for analyzing structural stress by utilizing stress 
resultants directly available from finite element codes or another source of similar data; 
and 

Figs. 1 1 and 12 illustrate experimental techniques for measuring structural 

stress. 
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DETAILED DESCRIPTION 
The following detailed description is presented according to the following subject 
matter outline, as delineated in the body of the description: 

Structural Stress Analysis by Enforcing Equilibrium Conditions and solving for a B , 
and a M 

Stress Analysis by Using Stress Distributions a x (y) and r^y). 
Stress Analysis by Using Stress Resultants. 
Special Applications - Partial Thickness Fatigue Crack. 
Special Applications - Non-Monotonic Through-Thickness Distributions. 

Calculation of Structural Stress by Conversion of Nodal Forces and Moments. 
Conversion of Nodal Forces and Moments Retrieved Directly from Shell Model. 
Conversion of Nodal Forces and Moments by Generating Stiffness Matrices and 
Nodal Displacements from the Shell Model. 

Conversion of Nodal Forces and Moments from Three-Dimensional Solid Model. 

Experimental Techniques for Measuring Structural Stress 
Monotonic Through-Thickness Distributions 
Non-monotonic Through-Thickness Stress Distributions 

Referring initially to Figs. 1A, 1B, 1C and 1D, it is noted that structures containing 
geometric discontinuities generally possess localized stress concentrations and 
corresponding through-thickness stress distributions. For example, the structure or 
25 plate 10 illustrated in Fig. 1 A includes a fillet weld 12 and defines a particular through- 
thickness stress distribution <r x (y). The stress distribution o x (y) is defined along a 
cross section of the structure 10 in a fatigue-prone region or weld toe 14 of the weld 12 
and may be obtained, for example, from a finite element model of the structure 10. Fig. 
1B illustrates the corresponding structural stress definition in terms of two components 
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a M , a B of the structural stress a s in the localized fatigue-prone region 14. The 
structural stress definition illustrated in Fig. 1B is effectively a simple structural stress 
distribution in the form of membrane and bending components a M , a B that are 
equilibrium-equivalent to the local stress distributions in Fig. 1 A. 

Through-thickness stress distributions a x (y) with respect to potential failure 
planes for a structure 10' including lap fillet joints 12' are illustrated in Fig. 1C. In such 
joint types, depending on actual geometry and loading conditions, three possible 
fatigue-prone regions or failure planes exist: (1) weld root failure along the weld throat; 
(2) weld toe failure; and (3) weld root/toe failure. The through-thickness stress 
distributions 0 x (y) correspond to cross sections of the structure 10' taken along the 
different fatigue-prone regions of the structure 10'. The corresponding equilibrium- 
equivalent simple structural stress definitions are illustrated in Fig. 1D in terms of the 
two components a M , a B of the structural stress a s in the different fatigue-prone regions. 

Referring now to Figs. 2A and 2B, local normal stress distributions <r x (y) along a 
typical failure plane are illustrated for structures 10" including a spot weld 12" (e.g., a 
resistance spot weld) or a laser lap weld 1 2". Regardless of whether the structure 1 0" 
includes a laser lap weld or resistance spot weld, the local normal stresses and 
structural stresses are essentially defined the same. Note that for resistance spot 
welds, see Fig. 2A, both the local stresses and the structural stresses are a function of 
circumferential angle. In the case of typical laser lap welds, see Fig. 2B. In the case of 
a laser lap weld, the weld area will occupy the entire cross-section without the half-width 
(w/2). If the stress quantities can be adequately defined in the x-y plane of Fig. 2B 
despite the fact that the long lap joint extends in a direction orthogonal to the x-y plane. 

Figs. 3 illustrates in further detail the local normal stress distribution a x (y) and 

the shear stress distribution r^fy) of a structure 20 or, more specifically, a structural 
member 20, in a fatigue-prone region 24 in the vicinity of a weld 22 of the structure 20. 
The values of o x (y) and r xy (y) are taken along line B-B of Fig. 3 for calculation 
purposes only. The value of t corresponds to the thickness of the structure 20 in the 
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selected cross section. The structural mid-plane 26 is illustrated in Figs. 3 and 4. Fig. 4 
represents a structural stress definition corresponding to a x (y) and r xy (y) along section 
A-A of Fig. 3 and includes, more particularly, the two components a M , a B of the 
structural stress a s in the localized fatigue-prone region 24. Also illustrated in Figs. 3 
5 and 4 is the variable d which represents the finite element size of a finite element model 
utilized in calculating structural stress according to the present invention. Typically, a 
single element row corresponds to the distance between sections A-A and B-B. It is 
noted that a B assumes a value of zero at the structural mid-plane 26. 

10 ^5 1 ■ Structural Stress Analysis by Enforcing Equilibrium Conditions and solving for o B . 
and a M 

¥U 

U1 A finite element model of a structure can be used to retrieve either stress outputs 

pj or nodal forces and displacements of the structure. If nodal forces and displacements 
f . are retrieved from the model, they can be used to calculate structural stresses in an 

IssSi 

3 

15 H- accurate manner. In using two and three-dimensional solid element models, it can be 
more convenient to retrieve stress outputs directly from the finite element model. 
Q Stress outputs retrieved in this manner must, however, be processed by enforcing 
equilibrium conditions at a reference location in the structure. 

According to the methodology of one embodiment of the present invention, 
20 structural stress a s in a localized fatigue-prone region of a structure (e.g., section A-A of 
Figs. 3 and 4) is defined in terms of the bending and membrane components <t b , a M oi 
the structural stress <r s in the region. The structural stress is calculated by enforcing 
equilibrium conditions in the region and performing a selected integration utilizing stress 
distributions a x (y) and r^fy) derived from an established source, such as a finite 
25 element model. 

As will be appreciated by those practicing the present invention, finite element 
models of a variety of structures have been, and may be, generated through computer- 
aided engineering. These finite element models may be used to derive individual 
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stresses or internal loads for an area of interest within the structure. Accordingly, the 
stress distributions ajy) and r xy (y) are typically derived from a finite element model 
and are described herein in the context of finite element models. It is important to note, 
however, that the stress distributions may be derived or input from any of a variety of 
sources, either existing or yet to be developed. For example, a number of finite 
difference models and so-called meshless calculation systems are currently under 
development for structural stress modeling. 

The structural stress distribution <r x (y) must satisfy equilibrium conditions within 
the context of elementary structural mechanics theory at both the hypothetical crack 
plane (e.g., at weld toe 14 in Fig. 1A) and a nearby reference plane, on which local 
stress distributions are known a priori from typical finite element solutions. The 
uniqueness of such a structural stress solution can be argued by considering the fact 
that the compatibility conditions of the corresponding finite element solutions are 
maintained at this location in such a calculation. 

It should be noted that in typical finite element based stress analysis, the stress 
values within some distance from the weld toe can change significantly as the finite 
element mesh design changes. This behavior is referred to herein as mesh-size 
sensitivity. While local stresses near a notch are mesh-size sensitive due to the 
asymptotic singularity behavior as a notch position is approached, the imposition of the 
equilibrium conditions in the context of elementary structural mechanics according to 
the present invention eliminates or minimizes the mesh-size sensitivity in structural 
stress calculations. This is due to the fact that the local stress concentration close to a 
notch is dominated by self-equilibrating stress distribution. 

1.1. Stress Analysis by Using Stress Distributions a x (y) and r v . (y). 

According to one method of the present invention, a cross section of interest 
(e.g., section B-B of Figs. 3 and 4), the location of which is defined by a unit width, is 
selected such that the localized, fatigue-prone, region lies adjacent to the selected 
cross section. For the purposes of describing and defining the present invention, it is 
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noted that elements that are separated by a single element distance from a selected 
cross section are said to lie adjacent to the cross section. 

A first component a M of the structural stress a s in the localized region (e.g., 
section A-A of Figs. 3 and 4) represents what may be referred to as a membrane 
component of the structural stress a s . A second component a B of the structural stress 
a s in the localized region represents what may be referred to as a bending component 
of the structural stress a s . The structural stress a s is calculated by combining the first 
component <r M of the structural stress and the second component a B of the structural 
stress as follows: 



By imposing equilibrium conditions between Sections A-A and B-B, the structural 
stress components a m and a b may be determined according to the following equations. 
Specifically, the first component <7 M of the structural stress a s in the localized region is 
determined according to the following equation: 



The stress distribution a x (y) along the selected cross section represents the local 
structural through-thickness stress distribution and is determined from the finite element 
model of the structure. The second component a B of the structural stress a s in the 
localized region is determined by solving the following equation for a B : 



B ' 




0 




0 0 



t t 
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The term y corresponds to a distance along the y-axis from y=0 to a material point of 
interest in the selected cross section. The term t corresponds to the thickness of the 
structure in the selected cross section and S represents the element size of the finite 
element model. The function r xy (y) represents the shear stress distribution and is 
5 determined from the finite element model of the structure. It is contemplated that, in 
practice, a row of elements may be used between the localized fatigue-prone region of 
a structure (e.g., section A-A of Figs. 3 and 4) and the selected cross section (e.g., 
section B-B of Figs. 3 and 4). 



10 ip the force balances in the x direction, evaluated along B-B. The equation for the second 
jjj component a B of the structural stress a s represents the moment balances with respect 
^ to Section A-A at y = 0. The integral term on the right hand side of the equation for the 



W? second component a B represents the transverse shear force as a component of the 

I structural stress definition. Accordingly, it is noted that, under certain conditions, such 

15 [T as if the transverse shear effects are not significant or if the finite element formulation 

IF" 

N provides adequate consideration of shear locking effects (particularly when using 

y * 

q relatively large size elements), the structural stresses can be calculated directly from 

^ section A-A as follows: 



The equation for the first component a M of the structural stress a s represents 



20 




0 
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1-2. Stress Analysis bv Using Stress Resultants. 

Referring briefly to Fig, 10, stresses and nodal quantities from shell or plate finite 
element models are often defined in a global coordinate system (x, y, z), depending on 
the finite element codes used. Given the definition of the structural stress components 
5 in Eq. 1 , it is the local coordinate system (x',y',z') that is convenient for calculating the 
structural stresses with respect to a weld, with local x' and y' being perpendicular and 
parallel to the weld direction, respectively. Consistent with the solid element model 
approach (e.g., see Fig. 3), three components of the stress resultants (sectional forces 
and moments), i.e., f x ,f z >, and m y , , at Section B-B in Fig. 10 can be used to calculate 
10 y3 the structural stress components at Section A-A: it is noted that structural stress may 
j| also be analyzed by utilizing stress resultants directly available from finite element 

fU codes or another source of similar data. Specifically, three components of the stress 

Ul 

Ul resultants f^f^ and m y . in at Section B-B in Fig. 10 can be used to calculate the 

fU 

_ structural stress components at Section A-A according to the following equation: 

^ 6(w y + 8f z .) 
15 m ° s = Gm + ° B = ~T + T 2 

where 3 and t represent the dimensional values illustrated in Fig. 10 and f x < t f^ and m y , 
represent the stress resultants illustrated in Fig. 10. 

According to this aspect of the present invention, a finite element formulation 
with six degrees of freedom at each node is assumed, i.e., six components of 
20 generalized forces at each node (three translational and three rotational). If stresses in 
the global coordinate system are used, they must be transformed to the local 
coordinate system before structural stress is calculated in the manner described above. 



25 1.3. Special Applications - Partial Thickness Fatigue Crack. 

Referring to Figs. 5A, 5B, and 5C, an alternative scheme is described herein for 
calculating structural stress a s where the selected cross section of the structure at issue 
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includes a partial thickness fatigue crack extending a distance tj below a surface of the 
structure. According to this alternative scheme, the stress calculation scheme 
discussed above with respect to Figs. 3 and 4, is altered by introducing two additional 
structural stress parameters a m f and c b * to enable the calculation of structural stress a s 
5 where the selected cross section of the structure at issue includes the partial thickness 
fatigue crack extending a distance tj below a surface of the structure 20. The 
parameters a m ' and a b f are illustrated in Fig. 5C. The overall through-thickness 
structural stress components a m and a b , illustrated in Fig. 5B, are calculated as 
^ discussed above with respect to the components a M and a B of Figs. 3 and 4. For the 
1 0 /? scheme of Figs. 5A, 5B, and 5C, the structural stress a s is calculated by combining a 



redefined first component a M of the structural stress and a redefined second 
component a B of the structural stress as follows: 



"™ According to this alternative scheme, the first and second components a M , a B of the 
15 CP structural stress a s in the localized region, with respect to t 1 are determined by solving a 
jS first equation through direct integration and by solving simultaneously three equations 
with three unknowns. 

A first equation of the four equations defines the sub-component a m 9 of the 
structural stress a s . A second equation of the four equations defines a force equilibrium 
20 condition. A third equation of the four equations defines a moment equilibrium 

condition. A fourth equation of the four equations defines a stress continuity condition. 

Regarding the first equation, since a m r is not mesh-size sensitive, it can be 
calculated by direct integration along B-B within the length of t-tj as follows: 



25 
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1 l ' h 1 5 

1 M 0 [ MO 



where <r m ' and a b ' comprise respective sub-components of the structural stress <r s and 
are taken relative to a reference line 26' of the structure 20. The position of the 
reference line 26' corresponds to the centroid or neutral axis of the region defined 
between t=0 and t=t~t } . The structural mid-plane 26 is also illustrated. 

By enforcing force balances along x and y, and moment balance at y=0 at 
n Section A-A, the second, third, and fourth equations are available. The second 
2 equation, defining the force equilibrium condition along the x-axis, is as follows: 



3 = 

If" 



1 I s 
where ° m = ~ \° x (y)<ty + ~ K*0>& . 

MO MO 



The third equation, defining the moment equilibrium condition is as follows: 



{°H-*B)t{t-\ } 



+ <x B t x 
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where a b is determined by solving the following equation for a b : 
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where the term t corresponds to the thickness of the structure in the selected cross 
section and 3 represents the element size of the finite element model. The function r xy 
(y) represents the shear stress distribution and is determined from the finite element 
model of the structure. 

5 The fourth equation, defining the stress continuity condition at y - t - t x is as 

follows: 

^ These three additional equations, defining three unknowns (a Mt <r Bt <r b r ), may be 

P 

*D solved simultaneously, according to conventional methods, to determine the values of 

10 i0 the first and second components o M , a B of the structural stress a s in the localized 

m 

m region. 

m 
ru 

3 

y, 1.4. Special Applications - Non-Monotonic Through-Thickness Distributions. 

15 Referring now to Figs. 6A, 6B, and 6C, a process for calculating structural stress 

y a 

O in a localized fatigue-prone region of a structure 30 defining a non-monotonic through 
^ thickness stress distribution is illustrated. The structure 30 includes a pair of welds 32 
defining weld toes 34 and is merely illustrated partially and schematically in Figs. 6A, 
6B, and 6C. Fig. 6A illustrates local normal and shear stress distributions of the 
20 structure 30 at the weld toe 34 of the structure 30. Fig. 6B illustrates respective local 
normal a x (y) and shear T xy (y) stress distributions at a reference section of the 
structure 30. Fig. 6C illustrates structural stress components relative to both upper and 
lower portions of the structure 30. 

The stress distribution c x (y) of Figs. 6A and 6B defines a minimum axial stress 
25 along a secondary axis 36* of the structure 30 where the transverse shear stress 
distribution r^fy) changes signs. The secondary axis 36' of the structure lies a 
distance t 2 below a surface of the structure 30. The structural mid-plane 36 is also 
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illustrated in Figs. 6A, 6B, and 6C. The cross-section references A-A and B-B are 
similar to those described above with reference to Figs. 3 and 4. The parameter t 2 can 
be determined based on the position at which the transverse shear stress changes 
direction, if there is no specified crack depth as a failure criterion. 

Two additional structural stress parameters a m f and a b f are illustrated in Fig. 6C. 
Fig. 6C illustrates these further sub-components a m ' and o b * of structural stress a s 
taken relative to the secondary axis 36* of the structure 30. Fig. 6C also illustrates the 
two primary components a M and a B of the structural stress a s where 

°s = ° M + G B 



10 ^0 and where the stress continuity condition at y = t - t 2 is as follows: 



Fig. 6C also illustrates the structural stress definitions corresponding to the non- 
p monotonic through-thickness stress distribution illustrated in Fig. 6A. If the stress 
1 5^" distribution a x (y) is symmetric with respect to the mid-plane 36 of the structural 

member 30, then t 2 = t/2 and the reference plane 36' coincides with the mid-plane 36. 
Generally, a M , a B , a b \ and a m f can be solved in a manner similar to the scheme 
illustrated above with respect to the embodiment of Figs. 5A, 5B, and 5C, provided t 2 
has been determined. The parameter t 2 can be determined either by considering 
20 equilibrium conditions with respect to the entire through-thickness section (t) or, more 
conveniently, by identifying the point at which the transverse shear stress distribution 
T xy(y) changes signs, as is also discussed above with reference to the location of the 
structural mid-plane 26 of Fig. 5B. 

Where the top surface of the structure 30 corresponds to the location of the 
25 maximum structural stress, the structural stress a s in the localized region may be 
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calculated by solving simultaneously the following equations for the two unknowns, o M , 



'2 t-h 



tter 



l 2 



°M + 



l 2 



t-h 



5 jy The two additional structural stress parameters a m f and a b f can be calculated according 
to the following equations: 



O 



1 r 

° m =t—t \° x (y)<ty 

1 l 2 0 



10 ff -^^ +(r »^^ = l^W^ 



For the calculation of structural stress for resistance spot and laser lap welds (see Figs. 
2A and 2B), if two or three-dimensional solid element models are used, the calculation 
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procedures are similar to the procedures described above with reference to Fig. 4. The 
fatigue crack typically initiates at the weld edge at the interface between the two sheets, 
where structural stress peak is located. A special case arises where the joint 
configuration and loading are symmetric with respect the neutral axis. In this case, the 
5 shear stress on the cross section along the symmetry line is zero and structural stress 
components may be calculated by substituting t 2 - 1/2 and r yx - 0. 



2.0. Calculation of Structural Stress bv Conversion of Nodal Forces and Moments. 
10 □ According to another embodiment of the present invention, structural stress a s in 

dg a localized fatigue-prone region of a structure is also calculated. This embodiment of 

^ the present invention involves conversion of relevant nodal forces and moments, or 

Wx internal force and moments at the nodes of a group of elements aligned with a weld toe 

m line, to sectional forces and moments and may or may not involve a direct integration 

15", step similar to that described herein with reference to Figs. 3-6. 

h As will be appreciated by those practicing the present invention, finite element 

u 

m models of a variety of structures have been, and may be, generated through computer- 
™ aided engineering. These finite element models may be used to derive individual 
stresses for an area of interest within the structure. Accordingly, the nodal forces and 
20 moments are typically derived from a finite element model and are described herein in 
the context of finite element models. It is important to note, however, that the nodal 
forces and moments may be derived or input from any of a variety of sources, either 
existing or yet to be developed. For example, a number of finite difference models and 
so-called meshless calculation systems are currently under development for structural 
25 modeling. 

2.1. Conversion of Nodal Forces and Moments Retrieved Directly from Shell Model. 
Specifically, referring to Figs. 7 and 8, according to one aspect of the present 

invention, structural stress a s in a localized fatigue-prone region of a structure is 
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calculated from a finite element shell model of the structure. Typically, as is illustrated 
in Fig. 8, the structure and the welds are modeled with four-node (quadrilateral) shell or 
plate elements 42. Initially, a shell element model 40, see Fig. 8, of the structure at 
issue is generated utilizing a conventional computer aided engineering application (see 
5 step 100). Finite element analysis is performed on the shell element model (see step 
102) and a determination is made as to whether nodal forces and moments for local 
elements of interest 42 can be directly retrieved from the output of the finite element 
analysis (see step 104). As is illustrated in Fig. 8, the elements 42 are positioned 
adjacent to a localized fatigue-prone region of a weld 44 and, as such, are identified for 
10 ^ structural stress extraction. If the nodal force and moment vectors for the local 
*0 elements 42 may be retrieved directly from the finite element shell model 40 then 

aft 

.J selected ones of the nodal force and moment vectors are converted to sectional force 

^ vectors n and moment vectors m with an appropriate mapping function that provides 

m 

W the consistent nodal loads, or work equivalent (see step 112). The conversion is 
1 5 ~z performed in a work equivalent manner with respect to said nodal displacements 

determined for said nodal force and moment vectors (see step 113). Specifically, the 
M> conversion is performed such that a quantity of work corresponding to said nodal 

en 

h displacements and nodal force and moment vectors is equivalent to a quantity of work 
^ for said nodal displacements and said sectional force and moment vectors n and m. 

20 The mapping function ideally should be the same as that used in the finite 

element program used in performing the structural analysis. The mapping function 
provides sectional forces (force per unit length) and moments (moment per unit length) 
corresponding to the nodal forces and moments. If high stress gradients exist, after 
mapping, the sectional forces and moments for a group of elements along the weld may 

25 not be continuous from one element to another. However, the total nodal forces and 
moments (as well as displacements) shared by neighboring elements at a nodal point 
or a grid point are correct from typical finite element solutions. The continuity in 
sectional forces and moments must be enforced by considering the neighboring 
elements and solving a system of linear equations formulated by using an appropriate 

30 mapping function for each element. The unknowns in the system of linear equations 
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are sectional forces and moments at the nodes adjoining neighboring elements. The 
system of linear equations are formulated using nodal forces and moments (outputs 
from a finite element shell or plate models) and element geometric information. After 
solving the system of linear equations, the sectional forces and moments in the 
direction of interest can be used to calculate the structural stresses to be described 
below. If stress gradients along weld lines are not high and element shapes are not 
distorted, the structural stresses can be directly calculated on an e!ement-by-element 
basis without resort to solving a system of linear equations, as discussed below with 
reference to the equation for structural stress a s . 

Structural stress c s in a localized fatigue-prone region of the structure may then 
be calculated according to the equation 

where a B is proportional to the sectional moment vector m and a m is proportional to the 
sectional force vector n. More specifically, the structural stress o s may be calculated 
utilizing the following equation: 

Ylmz n 
a s = a B + a M = — — + y 

where t corresponds to the thickness of the structure in the fatigue-prone region and z 
ranges from +t/2 at the top surface of the structure to -t/2 at the bottom surface of the 
structure. 

2.2. Conversion of Nodal Forces and Moments Using Stiffness Matrices and Nodal 
Displacements from the Shell Model. 

In some applications, the reference section B-B in Fig. 10 may not be easily 
defined. This situation arises if welds are rather close to each other or load transfer at 
a weld of interest is very localized. If the element sectional forces and moments (with 
respect to the reference element in Fig. 10) at Section A-A are available from a finite 
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element solution, the equilibrium requirements are automatically satisfied within the 
accuracy of the finite element solutions. In this context, if the nodal force and moment 
vectors can not be retrieved directly from the finite element model of the structure (see 
step 104) then the nodal force and moment vectors are computed by generating 
5 stiffness matrices and nodal displacements for the local elements in the fatigue-prone 
region from the finite element model (see step 106). The nodal force and moment 
vectors for elements of interest may then be computed by multiplying the stiffness 
matrices and nodal displacements to obtain global nodal force and moment vectors at 
nodal points of said local elements and transforming the resulting global force and 
1 0 p moment vectors from the global coordinate system to the local coordinate system of an 
element of interest (see step 1 08). As is noted above, the transformation is performed 

gg in a work equivalent manner with respect to said nodal displacements (see step 109). 

f~ More specifically, the global stiffness matrix and nodal displacements may be 

In used to determine local nodal force and moment vectors of an element of interest 
1 5 „ because the global element stiffness matrices and nodal displacements of the structure 

f7 at issue are known (see step 110). The structural stress is calculated as follows: 

M> 

S f 6(m y ) 

where t represents the dimensional value illustrated in Fig. 10 and f x , and m y , 

20 represent the stress resultants illustrated in Fig. 10. 

It should be noted that, according to this aspect of the present invention, the 
finite element software used must allow for generation of the stiffness matrices and the 
nodal displacements for the elements of interest. It should also be noted that, since the 
element nodal forces and moments are based on finite element solutions, the 

25 transverse shear effects considered earlier in the context of two-dimensional solid 
models are already taken into account in the shell element approach. 
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2.3. Conversion of Nodal Forces and Moments from Three-Dimensional Solid Model. 

As is illustrated in Fig. 9, structural stress a s in a localized fatigue-prone region of 
a structure may also be calculated from a finite element solid model of the structure by 
using noda! forces in a manner similar to that described above in the context of the 
shell/plate element models (see Figs. 7 and 8). A general three dimensional finite 
element solid model 40' including local elements 42' and a weld 44 is illustrated in Fig. 
9. A fatigue-prone plane of interest is illustrated along section A-B-C. Nodal forces on 
element faces can be either determined directly or through the local element stiffness 
matrix in a manner similar to that described in the context of Fig. 7. Once the nodal 
forces at the nodal positions are obtained, equivalent sectional forces and moments 
along Section A-B-C can then be obtained using consistent mapping or shape 
functions. The structural stress can then calculated with the equivalent sectional forces 
and moments in the manner described with reference to Figs. 7 and 8. In addition, 
equivalent transverse shear effects need to be considered as discussed earlier with 
respect to the two-dimensional solid models. 

It is contemplated that the various calculation schemes of the present invention 
are best performed with the aid of a computer programmed to execute the variety of 
calculation steps described herein or a computer arranged to read a suitable program 
stored on an appropriate recording medium. 

3. Experimental Techniques for Measuring Structural Stress 

Referring now to Figs. 1 1 and 12, it is noted that the structural stress calculation 
schemes of the present invention may also be utilized in fashioning experimental 
techniques for measuring structural stress by enforcing equilibrium conditions within a 
strip element adjacent to the weld toe of interest. 

3.1. Monotonic Through-Thickness Distributions 

Fig. 1 1 illustrates the local normal stress distribution a x (y) of a structure 20 or, 
more specifically, a structural member 20, in a fatigue-prone region 24 in the vicinity of 
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a weld 22 of the structure 20. The value of ajy) is taken along line A-A. The value of 
t corresponds to the thickness of the structure 20 in the selected cross section. The 
structural mid-plane 26 is illustrated in Fig. 1 1 . Strain gauges 60, or other displacement 
measuring devices, are arranged along cross-sectional lines B-B and C-C. 

Where through-thickness normal stress distributions take the form as shown in 
Fig. 11, stresses at the strain gauges can be calculated using strain gauge readings 
from measured positions along sections B-B and C-C. If the through-thickness stress 
distributions at the measured positions are approximately linear then the stress 
measurements can be decomposed as follows: 



where the superscripts B and C indicate the sections at which the stresses are 
measured and subscripts T and B refer the values at the top and bottom surfaces, 
respectively. Linearity of the stress distributions at the measure positions can be 
readily confirmed by using one or more strain gauge readings, either at the top or 
bottom, to see if a linear distribution along the surface is developed. 

The transverse shear resultant force effects at B-B, as described earlier in the 
context of finite element models, can be approximated by 



where / is the moment of inertia of the structural member 20. If the distance between 
sections A-A and C-C is small, the increase in the bending stress component at section 
A-A can be approximated as: 
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°b = °Z + 



L 



I 



The structural stress, as defined previously in Fig. 4, can then be approximated based 
on the measurements at sections B-B and C-C as: 



The distances among sections A-A, B-B, and C-C are typically measured in terms of 
fractions of thickness t. 

3.2. Non-monotonic Throuah-Thickness Stress Distributions 

Referring now to Fig. 12, it is noted that non-monotonic stress distributions can 
be handled in a manner similar to that described above with reference to Fig. 1 1 . Fig. 
12 illustrates the local normal stress distribution <r x (y) of a structure 30 or, more 
specifically, a structural member 30, in a fatigue-prone region 44 in the vicinity of a weld 
32 of the structure 30. The value of c x (y) is taken along line A-A. The value of / 
corresponds to the thickness of the structure 30 in the selected cross section. The 
structural mid-plane 36 is illustrated in Fig. 1 1 . Strain gauges 60 are arranged along 
cross-sectional lines B-B, C-C, and D-D. 

Measurements at section D-D, where the stress distribution <r x (y) 9 becomes 
approximately linear, provide the through-thickness mean stress a m as follows: 



Typically, a m =F/A (load F and cross-sectional area A are typically available). 
Accordingly, the structural stress component a M , as defined in Fig. 6C, becomes g^g, 
and 





2 
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B B 
G b - G r ~ (J 

G b — G T — G 



m 



m 



The structural stresses with respect to the upper half of the thickness (as shown in Fig. 6C) 
can be approximated as: 



It is noted that general monotonic and symmetric stress distributions can be 
handled in a manner to that illustrated in Figs. 1 1 and 12, with the exception that t 2 and 
the corresponding appropriate equilibrium conditions must be accounted for in the 
manner discussed above with respect to Fig. 6C. A structure including a partial 
thickness fatigue crack, as illustrated with reference to Figs. 5A-5C above, can also be 
treated with a technique similar to that illustrated with reference to Fig. 1 1 , as long as 
the two strain gauges 60 utilized at the bottom of the structure 20 (see Fig. 1 1 ) are 
positioned a distance t, from the top edge (see Fig. 5B). The specific calculation 
procedures are then the same as those illustrated with reference to the monotonic 
through-thickness distribution case of Fig. 11. 

It is also noted that the present application discusses conventional strain gauges 
as a suitable means for measuring displacement/strain merely for the purposes of 
illustrating the present invention. It is contemplated that a variety of displacement/strain 
measurement devices may be employed according to the present invention including, 
but not limited to, imaging devices, laser-based devices, fiber-optic gauges, ultrasonic 
gauges, etc. 

For the purposes of describing and defining the present invention it is noted that 
the term "substantially" is utilized herein to represent the inherent degree of uncertainty 
that may be attributed to any quantitative comparison, value, measurement, or other 
representation. The term "substantially" is also utilized herein to represent the degree 
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by which a quantitative representation may vary from a stated reference without 
resulting in a change in the basic function of the subject matter at issue. 

As will be appreciated by those practicing the present invention, it is possible to 
utilize alternatives to the selected integrations described in detail herein. Indeed, any 
5 mathematical or procedural operation that has a result that is substantially equivalent to 
the result of a selected integration may be performed according to the present 
invention. For the purpose of defining the present invention, a claim recitation of an 
operation having a result substantially equivalent to a result of a given integration reads 
on the integration itself and the equivalent operation. 
10^ Having described the invention in detail and by reference to preferred 

>D embodiments thereof, it will be apparent that modifications and variations are possible 

without departing from the scope of the invention defined in the appended claims. More 
j^j specifically, although some aspects of the present invention are identified herein as 
tq preferred or particularly advantageous, it is contemplated that the present invention is 



15^" not necessarily limited to these preferred aspects of the invention. 
f* 4 What is claimed is: 



- 28 of 51 - 



